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ABSTRACT 

Wc estimate the contribution to the extragalactic gamma-ray background (EGRB) 
from both intrinsic and cascade emissions produced by blazars using a simple semi- 
analysis method for two models of the blazar gamma-ray luminosity function (GLF). 
For the cascade emission, we consider two possible contributions: one is due to that 
the flux of the cascade emission is lower than the Fermi LAT sensitivity (case I), which 
is independent on the extragalactic magnetic field (EGMF), another is due to the fact 
that the flux of the cascade emission is larger than the Fermi LAT sensitivity but the 
emission angle is larger than LAT point-spread- function (PSF) angle (case II), which 
depends on the EGMF. Our results indicate that (1) blazar contribution to the EGRB 
is dominant although it depends on the GLF model and the EGMF; (2) the EGMF 
plays an important role in estimating the contribution from the cascade emission 
produced by blazars, the contribution from the cascade emission will significantly 
alter the EGRB spectrum when the strength of the EGMF is large enough ( say 
-Begmf > 10^^^ G); and (3) since the cascade emission in case II reaches a saturation 
when the strength of the EGMF is ^ 10~^^ G, it is very possible that the contribution 
from the cascade emission produced by blazars can be considered as another method 
to probe the upper limit of the strength of the EGMF. 
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1 INTRODUCTION 

The gamma-ray sky observed by the Large Area Telescope 
(LAT) on board the Fermi Gamma Ray Space Telescope 
consists of the resolved gamma-ray emitters, such as nor- 
mal galaxies, active galactic nuclei (AGNs), gamma-ray 
bursts (GRBs), and pulsars etc., and the diffuse gamma- 
ray radiation including emission from the Galaxy and 
the extragalactic gamma-ray background (EGRB) (e.g., 
lAbdo et al'l l2010al fbl). The origin of the EGRB is poor 
known, which is one of the fundamental unsolved prob- 
lems in astrophysics. The EGRB was first detected by 
the SAS-2 satellite (|Fichtel et all Il975l ). Latter, in the 
1990s the Energetic Gamma-Ray Experiment Telescope 
(EGRB) on board the Compton Observatory measuerd 
its spectrum at 0.03 - 50 GeV with a good accuracy 
jSreekumar et al.lll998l : IStrong et al.ll2004l ). Recently, LAT 
has made a new measurement of the EGRB spectrum 
lAbdo et all l2010bl ). This has been found to be consis- 
tent with a featureless power-law with a photon index 
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of ~2.41 in the 0.2 - 100 GeV energy range. The ob- 
served integrated flux above 100 MeV is (1.03 ± 0.17) x 
10~^ photons cm~^ s~^ sr~^, which is lower than the one 
of (1.14 ± 0.05) X 10~^ photons cm~^ s~^ sr~^ observed by 
EGRET (|Strong et al.ll2004l : lAbdo et al.ll2010bl ). 



Various kinds of 
have been proposed 
to t h e EGRB , such 
19931: IChiang 



unresolved gamma-ray sources 
as the possible_ contributors 

IS AGNs 

Mukheried 1 19981: 



Padovani et al.l 



Giommi et al.l 



Narumoto fc Totanill2006l: iDermeilWoTl: ICao fc Bai I 



Inoue fc Totanil l2009l: IVenters et all l2009l: IVenterd I2OI0I: 



2006 



2008 



Inouell201ll: lAie llo ot al.l 201 ll). star- fo rming galaxi es (e.g. 



Pavlidou fc Fields ,2003 : iFields et all .2010: Makiva et all 



201lf). and starbur st galaxies (e.g.. [Thompson et al.ll2007l : 



Makiva et al.|[201ll ). Since blazars are the dominant extra- 
galactic gamma-ray sources, it is naturally expected that 
the intrinsic emission of an unresolved population of blazars 
would account for a sizable contribution to the EGRB. 
However, because of the uncertainties of gamma-ray lumi- 
nosity function (GLF) and the spectrum index distribution 
(SID) of blazars, this con tribution to the EGRB from 
blazars has been debated. IStecker fc Venter^ (|201ll ) and 
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Ilnoue fc Totanil (|2009l ) suggested that this contribution is 
dominated by the emissi on from unr e solved blazars, while 
lAiello et all (j201lh and lAbdo et all (j2010j l argued that 
unresolved blazars produce a small fraction of the EGRB. 

Additionally, very high energy (VHE) photons from 
blazars will interact with the UV-IR photons of extragalac- 
tic background light (EBL) to produce pairs of electrons 
and positrons, and the pairs will boost cosmic microwave 
background (CMB) photons to gamma-ray energies through 
inverse Comptori scattering proce s s (e.g . , IPai et alj I2OO2I; 
Yang et al.l I2OO8I : iTavecchio et all I2OI0I : iNeronov fc Vovkl 



2OIOI ). Cascade contributions to the EGRB have been con- 



sidered by rnany a u thors fe.g.. [Strone fc Wplfendaiel 1 19731 : 



IStrong et al] 1 19741 : iKachelriess et al.l 120111 )" The cascade 



flux produced by VHE photons of blazars, which is lower 
than the LAT sensitivity, will make contribution to the 
EGRB. Besides, for a large extragalactic magnetic field 
(EGMF) strength, the cascade photons will spread over a 
large angle, so that LAT cannot identify this emission from 
a point source, which gives rise to another contribution to 
the EGRB. The latter contribution depends on the strength 
of the EGMF. In the case of high collective high-energy in- 
tensities of blazars or high EBL intensit ies, the cascade con- 
tribu t ion to the EGRB is s i gnifica nt (jCoppi fc AharonianI 
ll997h . lKneiske fc MannheimI (20081 ) have estimated the con- 
tribution radiation from the first generation pairs to the 
EGRB, and they h ave found that t his c ascade contribu- 
tion is sizable. But Ilnoue fc Totanil (|2009( ) suggested that 
this ca scade contribution is negligible. Very recently. I Venters! 
l|2010h has also estimated this cascade contribution by us- 
ing a Monte Carlo program, and found that the cascade 
radiation greatly enhances the contr ibution to the EGRB 
from b lazar s. It should be n o ted th at iKneiske fc Mannheiml 
lj2008l ) and Ilnoue fc Totanil (j2009l ) have ta ken neither the 
SID of blazars nor the EGMF into account. IVentersI (|2010l ') 
included the effect of the SID of blazars but ignored the 
effect of the EGMF on the cascade contribution. 

In this paper, we revisit the cascade contribution to the 
EGRB from blazars by assuming that the observed LAT 
GeV spectrum could be extrapolated into TeV energy range 
with a power-law spectrum and the SID, in particular, we 
study the impact of the EGMF on the spectrum of the cas- 
cade radiation. In §2, we present the method we used to 
calculate the contribution from blazars and present our re- 
sults. Finally we give a brief discussion and conclusions in 
§3. 



2 CONTRIBUTION OF THE EGRB FROM 
BLAZARS 

2.1 Blazar gamma ray luminosity function 

To estimate the contribution to the EGRB from unresolved 
blazars. The blazar GLF is required. However, due to the 
small sample size, it is difficult to build a GLF directly us- 
ing the current gamma-ray- loud blazar sample. On the other 
hand, luminosity functions (LFs) of blazars at other wave- 
lengths (lik e radio or X-ray) are wid e ly studied in previou s 
works (e.g.. lDunlop fc Peacocljll990l : iHasinger et al.ll2005h . 
Moreover, it is well believed that the gamma-ray emission of 
blazars would correlate with the emissions at lower energy 



bands. For instance, it was found that there is a good cor- 
relation bet ween the gamma-ray luminosity and radio lum i- 
nosi ty (e.g.. lAck crman n et all I2OIII : iGhirlanda et al.l I2OIII ) , 
and IChisellini et al, ()2oTol )" found a positive correlation be- 
tween the jet power and the luminosity of the accretion disc 
in some blazars. Hence, the GLF is typically built from the 
luminosity functions in other wavelengths by using the corre- 
lation between gamma-ray emission and lower energy band 
emissions. 

For the purpose of our analysis, we will adopt 
the best-fit mode ls of blazar GLF constructed by 
iNarumoto fc Totanil ()2006h : pure luminosity evolution 
(PLE) a nd luminosity-dependent de nsity evolution (LDDE) 
models. INarumoto fc Totanil (|2006l ) limited the model pa- 
rameters by using likelihood analysis of the observed redshift 
and gamma-ray flux distributions of the EGRET blazars 
and found that the LDDE model gives a better fit to the 
observed distributions than the PLE model. For the blazar 
GLF p^{L^, z) an d its parameters in ea c h mo del, we totally 
use them given by INarumoto fc Totanil (|2006l ). 



2.2 The intensity of the EGRB from intrinsic 
emission of blazars 

In order to estimate the EGRB intensity from the intrin- 
sic emission of blazars, we need the SID of Fermi-LAT re- 
solved blazars. In the clean sample of the First LAT AGN 
Catalog (ILAC), corresponding to 11 months of data col- 
lected in science operation mode, there are 523 blazars 
jAbdo et al . 20J,0.a,). Their energy spectra between 0.1 GeV 
to 100 GeV can be fitted with a power-law spectral with 
photon index F. Their photon index distribution is com- 
pa tible with a Gaussian distribution. However, as argued 
bv lAbdo et al.l ()2010ah . sources with hard spectra are more 
easily detected by LAT, therefore the observed photon in- 
dex distri bution tends to hard er than the intrinsic one. As 
shown in lAbdo et al.l ()2010ah (see their Fig.l), when the 
integrated flux -F^,>ioo above 100 MeV is larger than 7 x 
10~* photon cm~^ s~^, LAT detected all sources irrespective 
of their photon indices, fluxes, or positions in the specific sky. 
For the sources with -Ft.,>ioo ^ 7 x 10^* photon cm"'^ s~^, 
their spectral photon index distribution, which is the intrin- 
sic SID, is a lso compatible with a Gaussian distribution and 
is given by (|Abdo et al.ll2010al ) 



(1) 



dN 
dT 

with a mean of /i = 2.40 ± 0.02 and a dispersion of cr = 
0.24 ±0.02 . 

When the GLF and the intrinsic SID are given, 
the EGRB intensity, Jegrb(^)i in units of photons 
cm~^ s~^ sr" 
is given by 



dV 



dN 
If 



dz 



d^V 
dzdQ 



(2) 



dL^p^Fj{E, z, L-y, V) 



-t(E,z) 



(1 - a;(F^,>ioo)) 



where Fmin = 1.2 and Fmax ~ 3.0 are minimum and maxi- 
mum values of the photon index, Zmax = 5.0 is the maximum 
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Figure 1. The o ptical depths expected in the EBL model of 
iFinke et al] boioh as a function of the photon energy E for 
sources located at z =0.02, 0.2, 0.5, 1.0, 1.5, 2.0, 3.0, 4.5 from 
bottom to top. 



redshift, 



10 erg s and 



erg s 



are the minimum and maximum luminosities respectively; 
F1^[E, z, L^,T) is the intrinsic photon flux at energy E of 
a blazar with gamma-ray luminosity L-y and a power-law 
spectrum at redshift z and is given by (jVenters et al.ll2009l l 



F;(£;,2,L^,r) = 



R2 
^0 



100 MeV' 



(3) 



where E^ = 100 MeV • 1.6 x 10""* erg, t{E,z) is the op- 
tical depth of the EBL for the sources at redshift z emit- 
ting gamma - ray p hoton energy E, here we use the model of 
iFinke et al.l (|2010l ) to derive r(E, z); The optical depths ex- 
pected in this EBL model as a function of the photon energy 
E for sources at different redshifts are presented in Fig[T] 
i^(f7,>ioo) is the detection efficiency of LAT at the pho- 
ton flux F^,>ioo , which corresponds to the integrated flux 
above 100 MeV from the sources at redshift z with and 
r, and J;.>ioo = F.^ (100 MeV,z,L^,r) x 100 MeV/(r- 1) 
(| Venters et al.l [20091 ): rj is the normalization factor and r] = 



where A'obs is the number of the sample, and A^o: 



the expected number and is given by 



4tv 



dV 



dN 

dr 



dz 



dzdQ 



(4) 



dL^ p^uj{Fj ^yioo) 



Here, the detect ion efficiency of LAT reported by 
lAbdo et all (|2010al ') and A'obs = 523 are used. 



2.3 The intensity of the EGRB from cascade 
emission of blazars 

For the cascade emission of blazars, we follow the geome - 
try of the cascade process given by iDermer et al.l (|201lf ) . 
In this geometry, VHE photons emitted at an angle di by 
a blazar at distance d could convert into electron-positron 
pairs via photon-photon absorption with the photons of the 
EBL after travelling a mean distance A-,.^ — d/T(E\ruE, z). 



When the pairs are deflected by the EGMF by an angle 
Sdfl, these pairs could scatter CMB photons to gamma- 
ray energies, where Odd = At/j^l for At < Acoh, and 
fdfl = AT/rL(Acoh/AT)^''^ for At > Acoh, At is the cooling 
length for electrons with energy 7 due to Thomson scat- 
tering, tl is the electron/positron Lar mor radius, and Aco h 
is the coherence length of the EGMF jPermer et aLlbOllh . 
The redirected secondary gamma-ray photons would arrive 
at an angle 9 — Xj^Oda/d to the line of sight. In calcu- 
lations. At ~ 3meC^/4crTf/cMB(l + z)"'7, where (jt is the 
Thomson cross section, and rrie is the rest mass of electron, 
and UcMB = 4 X 10~^^ erg cm~^ is the CMB energy den- 
stiy at z — 0.0; tl = meC^'y /cBegmf , where Begmf is the 
EGMF of strength and e is the elementary charge. Here, we 
set the coherence length of the EGMF Acoh ~ 1 Mpc. 

There are two possible cases for the contributions of 
the cascade emission to the EGRB. In the first case (case 
I), the cascade emission could make contribution to the 
EGRB if the flux of the cascade emission is lower than 
the LAT sensitivity. In the second case, although the flux 
of the cascade emission is larger than the LAT sensitivity, 
the angle between of the redirected secondary gamma-ray 
photons and the line of sight is larger than LAT point- 
spread-function (PSF) angle, i.e. 9 > SpsF, the cascade emis- 
sion also will not be attributed to a point source by LAT, 
and then make contribution to the E GRB, where 0psF = 
i^(0.001S)-°-^**[l + (0.001£;/15)2]°-^'' (jTavlor et al.ll201lh . 
In the following, we will calculate the contributions of cases 
I and II, respectively. 

To calculate the cascade flux, the intrinsic TeV flux is 
necessary, which is poor known due to the uncertainties of 
physics properties of blazars and the EBL. An approach to 
derive the possible intrinsic TeV spectrum is to extrapolate 
the observed LAT GeV spectrum to TeV energy range by 
using a power-law spectrum. Here, we consider the extrap- 
olated TeV spectrum as the intrinsic TeV spectrum, which 
is written as 



Fvue'''"'{Ev-he, z, 1/7, F) = 



L-,{l + z)^ 



-(■ 



E^ 



(5) 



47rd2(z)£;2 nOOMeV 

Therefore, the total cascade flux of inverse Compton- 
scattering CMB photons produced by the flrst generation 
of electrons through photon-photon pair production is given 
as (|Dermer et al.li201ll : iHuan et allboill ) 



T-,cascadc / rp 

^1 



8l7r 



/; 



16A3 (1 + z)^Ucmb 

c?7 



(6) 



lfCMB(l + ^) 



78exp(3ec/472ecMB(l + z)) 



r 



5.11 X 10" 

106 



z,L^,r){l-e 



where the dimensionless energy Ec = £-10V(5.11-10"), and 
ecMB = 1.24x 10~^ in rricC^ units is the average CMB photon 
energy at z = 0.0, and Ac = 2.426 x 10"^" cm is the Compton 
length. For the maximum dimensionless energy of intrinsic 
TeV photons emax, we take tmax = 2.0 x 10*, corresponding 
to iSvHE = 100 TeV. The emission from the secondary gener- 
ation electrons or more th an secondary generation elec trons 
is negligible at GeV range (jKneiske fc Mannheimll2008l l. The 
energy losses of the electrons by synchrotron emission are 
negligible since the ratio between the cooling times of in- 
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verse Compton scattering CMB photons and synchrotron 
emission is very small, ~ 10~^ even when the upper limit 
value of the EGMF ~ 10"^ G is used. 

The EGRB intensity in case I, independent of the 
EGMF, is given by 



'Fegrb(-E') 



dT 



dz 



dzdO, 



(7) 



dL^ Fj^^'^'^'^° {E , z, L-y, r) 



/t / T-icascadc\\ 

(1 cj(F^,>ioo )) , 



where F^f^Jfocf 



2F^'''"''"^''{E,z,L.,,r)W0 MeV 
•(100 MeV/£;)"^•^ which is derived by assuming that 
the cascad e spectrum is a power-law spectrum with photon 
index 1.5 (|Strong et a"Llll974l : pTavecchio et ai] |2010''). 

The EGRB intensity in case II, which depends on the 
EGMF, is written as. 



dz 



dzdO, 



(8) 



dL^PiF^'^eyi^^^ {E, z,L^,r) 



— t(E,z) / j-,cascadc\ 

xe 'io(F^^yioo ) • 

In case II, to avoid the situation that the produced pair are 
isotropized by EGMF when the EGMF of strength is large 
enough, we set GdR < ■k jl. 



2.4 Total EGRB intensity of blazars and its 
contribution to the EGRB 

As mentioned above, there are two possible cases (cases I 
and II) for the contributions of the cascade emission to the 
EGRB. If we only consider case I for the cascade emission, 
i.e., the flux of the cascade emission is lower than the LAT 
sensitivity, the total EGRB intensity of blazars is given by 



-Fegrb(-E) — -Fegrb(-B) + -Fegrb(-^) 



(9) 



The cascade emission is independent on the EGMF in this 
case. In fact, the contribution of the cascade emission to 
the EGRB would include the contributions in both cases 
I and II (in the latter the flux of the cascade emission is 
larger than the LAT sensitivity, however the angle between 
the redirected secondary gamma-ray photons and the line 
of sight is larger than LAT PSF angle) . Therefore, the total 
EGRB intensity of blazars is given by 

Fiomi(E) + [i^EiRB(-B) + F^^G^s.{Ey\ ■ (10) 



-Fegrb (-E) 



It should be noted that the cascade emission depends on the 
EGMF in this case. 

At first we calculate the intrinsic emission, the cascade 
emission, and the total EGRB intensities for the LDDE 
model of the blazar GLF and the calculated results are 
shown in Fig[2l In the top panel of Fig[2] we consider the 
case I of the blazar cascade emission, i.e., the total intensity 
is given by Eq. (9) and the cascade emission is independent 
on the EGMF; thin solid, dotted, and thick solid lines repre- 
sent the intrinsic, cascade, and total intensities, respectively. 
We have found that the contribution of the blazar cascade 
emission in case I to the EGRB could not be negligible, and 



would enhance the contribution of the blazar emission to 
the EGRB at high energies. At lower energies, however, this 
contribution is negligible. 

We now consider the contributions of the cascade emis- 
sion in both cases I and II, i.e., the total intensity is given 
by Eq. (10). Since the cascade emission in case II depends 
on the EGMF, we calculate the intensities with different 
strengthes of the EGMF. The results are shown in the lower 
panel of Fig[2] It can be seen from the figure that the inten- 
sity of the blazar cascade emission in case II increases with 
increasing of the strength of the EGMF, and will reach a 
saturation at some value of the EGMF. From our calcula- 
tions, the contribution of the blazar cascade emission in case 
II reaches a saturation when the strength of the EGMF is 
Begmf ~ 10"^^ G; the contribution of the blazar cascade 
emission with Begmf < 10^^** G in case II is negligible, 
and this contribution becomes significant compared with the 
contribution of the blazar intrinsic emission at several GeV 
to several tens GeV with Begmf > 10~^^ G. If the strength 
of the EGMF satisfies that Begmf > 10"^^ G, the total 
contribution of blazars can account for the large fraction of 
the EGRB at several GeV to several tens GeV. 

In FigO we present the results by using the PLE model 
of the blazar GLF. It is found that the contributions of both 
the blazar intrinsic emission and the cascade emission in 
case I are significantly larger than those derived by using 
the LDDE model of the blazar GLF. The features of the 
contribution from the blazar cascade emission in case II is 
similar with those in Fig. (2] If the strength of the EGMF 
is larger ( e.g., Begmf > lO^^'^ G), the contribution of the 
cascade emission to the EGRB is also significant. Assuming 
the PLF model, the total contribution of blazars can explain 
the bulk of the EGRB, irrespective of the strength of the 
EGMF. Moreover, it seems that the model of the blazar 
GLF has an impact on the contributions to the EGRB for 
both the cascade emission in case I and the blazar intrinsic 



3 DISCUSSION AND CONCLUSIONS 

We have studied the impact of the cascade radiation on 
the contribution of blazars to the EGRB by using a sim- 
ple semi-analytic model. In particular, we take the effect of 
the EGMF on the cascade contribution from blazars into 
account. We have found that if the strength of the EGMF 
is large enough (Begmf > 10^ ^'^ G), the cascade contribu- 
tion could significantly alter the spectrum of the EGRB at 
high energies, while with small strength of the EGMF large 
enough (Begmf < lO"'^* G), the cascade contribution is 
small, but cannot be ignored (see Figs. 2 and 3). 

There is little information about the strength of the 
EGMF up to now. Measurements of Faraday rotation of po- 
larization of radio radiation from distant quasars place up- 
per li mits on the EGMF strengths at the level of ~ 10~® G 
(e.g.. lKronb"er3ll994l : lBlasi et al.lll999l '). On the other hand, 
the possible cascade emission from resolved blazars provides 
a method to probe the lower bands on the EGMF. The 



lower bounds on the EGMF at the level of 10 



-10" 



were obtained by u s ing this method ( e .g., Tavlor et al.l 201 ll : 



Dermer et al.ll201ll:lHuan et al.|[201]]: iTavecchio et al.l 



Neronov fc Vovkll2010l : iDolag et al.ll201lh . In this pap 
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Figure 2. Blazar contribution to the EGRB assuming the LDDE 
model of the blazar GLF. Top: the calculation of the blazar in- 
trinsic contribution (thin solid line) and the cascade contribu- 
tion of case I (dotted). The thick solid line is the sum of the 
two contributions. Bottom: comparison of the blazar intrinsic 
contribution (thin solid line) and the cascade contributions in 
case I (dotted line) and in case II (dot-dashed lino ). For the 
cascade contribution of case II, we calculate with Begmf = 



10" 



',io-",io- 



, 10" 



, 10" 



G (dot-dashe line: from thin to 



thick), respectively. The thick solid line is the sum of the blazar 
intrinsic contribution and the cascade contributions of cases I and 
II with BegmfIO"" G. 




have shown that the contribution of the blazar cascade emis- 
sion to the EGRB very possibly provides another method to 
set the upper bands of the EGMF. In order to do so, the 
prerequisites are that we have the well determined GLF and 
we can well estimate the contributions of other sources to 
the EGRB. We have shown that a possible upper bound of 
the EGRB is ~ 10~" G. 

In estimating the cascade contribution to the EGRB , 
there is an uncertainty in the EBL model. IVentersI (|2010l ) 
have found that the amount of cascade radiation is sensitive 
to the EBL model. Actually, o ne of the EBL mode ls IVenterd 
(|2010D adopted is hig h level (TStecker et al1l2006l) . and the 
other one is low level (|Gilmore et al.ll2009l ). Ho wever, recent 
study indicated that the EBL intensity is low (|Abdo et al.l 
l2010dl ). and the commonly used several EBL models are 
low levels and have sim ilar results (e.g.. iFinke et af]|2010l : 
iDommguez et ahl 1201 ih . It is therefore expected that the 
impact of uncertainty of the EBL model on the contribution 
of the cascade emission is slight. 

Based on the observed gamma-ray sample, various 
models of the blazar GLF have been c o nstructed (e.g., 
Stecker fc Venters! 120111: lAiello et al.ll201ll : llnoue fc Totanil 



2009t ). lAiello et alj (|201lf ) have found that these models 



at redshift zero are relative consistent, while the signifi- 
cant discrepancy between these models will appear at red- 
shift one. The LDDE model we used here expects smaller 
number of blazars compared with th ose expected in th e 
mo dels of IStecker fc VeritersI (|201ll 'l. lAiello et all (l201ll) 
andllnoue fc Totanil (|2009l '). The model of lStecker fc Venterd 
(|201lF is similar with the PLE model we used in this paper. 
Actually, the determination of the blazar GLF solely based 
on the observed gamma-ray sample is not certain, due to 
the sources confusion effect and the active featu r e of b lazars 
(| Venters fc Pavlidou.|[2oTl[ ). IStecker fc Venters! (120111 ) indi- 
cated that the sources confusion would lead to the under- 
estimate of the contribution to the EGRB below 1 GeV 
from blazars. The active flare can bring the intrinsically faint 
blazars to the bright blazars sample, which will flatten the 
faint-end slope of the observed blazar source counts. This 
flare effect also can cause underestimate of the contribution. 
In this sense, we conclude that we cannot rule out blazar 
dominance of the EGRB from our calculations. 



10" 
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Figure 3. Same as Fig|2]except for the PLE model of the blazar 
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